Coulomb Drag of Massless Fermions in Graphene 
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Using a novel structure, consisting of two, independently contacted graphene single layers sepa- 
rated by an ultra-thin dielectric, we experimentally measure the Coulomb drag of massless fermions 
in graphene. At temperatures higher than 50 K, the Coulomb drag follows a temperature and car- 
rier density dependence consistent with the Fermi liquid regime. As the temperature is reduced, 
the Coulomb drag exhibits giant fluctuations with an increasing amplitude, thanks to the interplay 
between coherent transport in the graphene layer and interaction between the two layers. 

PACS numbers: 73.43.-f, 71.35.-y, 73.22.Gk 



Bilayer systems formed by two layers of carriers in 
close proximity are a fascinating testground for electron 
physics. In particular, the prospect of electron-hole pair 
(indirect exciton) formation, and dipolar superfluidity 1 
has fueled the research of electron-hole bilayers in 
GaAs/AlGaAs heterostructures^. Graphene^ is a par- 
ticularly interesting material to explore interacting bilay- 
ers. The symmetric conduction and valence bands, and 
the large Fermi energy favor correlated electron states at 
elevated temperatures^. The zero energy band-gap al- 
lows a seamless transition between electrons and holes in 
each layer, and obviates the large inter-layer electric field 
required to simultaneously induce electrons and holes in 
GaAs bilayers^. Coulomb drag, a direct measurement of 
inter-layer electron-electron scattering^ can provide in- 
sight into the ground state of two— and one— electron 
systems, as well as correlated bilayer state s 11 ! 12 . Here we 
demonstrate a novel, independently contacted graphene 
bilayer, and investigate the Coulomb drag in this system. 

Two main ingredients render the realization of inde- 
pendently contacted graphene bilayers challenging. First, 
an ultra-thin yet highly insulating dielectric is required 
to separate the two layers. Second, a method to position 
another graphene layer on a pre-existing device with min- 
imum or no degradation is needed to create the second 
layer of the structure investigated here. The fabrication 
of our independently contacted graphene bilayers is de- 
scribed in Fig. 1. First, the bottom graphene layer is 
mechanically exfoliated onto a 280 nm thick Si02 dielec- 
tric, thermally grown on a highly doped Si substrate. 
E-beam lithography, metal lift-off, and etching are used 
to define a Hall bar on the bottom layer [Fig. 1(a)]. A 
7 nm thick AI2O3 is then deposited on the bottom layer 
using a 2 nm oxidized Al interfacial layer, followed by 5 
nm of Al 2 03 atomic layer deposition^. The second, top 
graphene layer is also mechanically exfoliated on a similar 
Si0 2 /Si substrate. A poly methyl metacrylate (PMMA) 
film is applied on the top layer and cured. Using an 
NaOH etch 14 , the PMMA film along with the graphene 
layer, and the alignment marks are detached from the 
host substrate, forming a free standing membrane. The 
membrane is placed face down on the substrate contain- 
ing the bottom graphene layer [Fig. 1(b)], and aligned 



with it. A Hall bar is subsequently defined on the top 
layer [Fig. 1(c)]. Ten back-gated, independently con- 
tacted graphene bilayers have been fabricated and inves- 
tigated in this study, all with similar results. We focus 
here on data collected from three samples, labeled 1, 2, 
and 3, with mobilities between 4,200-12,000 cm 2 /Vs for 
the bottom layer, 4,500-22,000 cm 2 /Vs for the top layer, 
and with inter-layer resistances of 1-20 Gf2. These struc- 
tures are markedly different from graphene bilayers exfo- 
liated from natural graphite, consisting of two graphene 
monolayers in Bernal stacking^. 

We now turn to the individual layer characterization. 
The layer resistivities (p) and Hall densities measured 
for sample 1 at a temperature T=4.2 K, as a function 
of back-gate bias (Vbg) are shown in Fig. 2(a) and 
2(b), respectively. The potential of the both layers is 
held at zero (ground) for all measurements presented in 
this study. The bottom layer dependence on the applied 
Vbg shows ambipolar conduction and a finite resistance 
at the charge neutrality (Dirac) point, consistent with 
the expected response of gated monolayer graphene 1 ^. 
More interestingly, the top layer resistivity also changes 
as a result of the applied Vbg- This observation indicates 
an incomplete screening of the gate-induced electric field 
by the bottom layer—, which is most pronounced in the 
vicinity of the charge neutrality point, consequence of the 
reduced density of states in graphene. As we show below, 
we can quantitatively explain the layer resistivities and 
densities dependence on Vbg- 

Figure 2(c) shows the band diagram of the graphene 
bilayer at Vbg=® V; for simplicity the gate Fermi en- 
ergy and the charge neutrality point in the two layers 
are assumed to be at the same energy. Once a finite 
Vbg is applied, finite charge densities are induced in both 
top (tit) and bottom (rig) layers [Fig. 2(d)]. The dif- 
ference between the gate and bottom layer Fermi lev- 
els is distributed partly across the Si02 dielectric, and 
partly on the Fermi energy of the bottom graphene layer: 
eV BG = e 2 (n B +n T )/C S i0 2 + E F (n B ) (1); E F (n) = 
TiVFy/itn is the graphene Fermi energy measured with 
respect to the charge neutrality point at a carrier den- 
sity n, e is the electron charge, uf=1.1x10 6 m/s is the 
graphene Fermi velocity, and Csio 2 denotes the Si02 di- 
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FIG. 1: (color online) Optical micrographs (top) and schematic representation (bottom) of the fabrication process of an 
independently contacted graphene bilayer. (a) Hall bar device fabrication on the bottom graphene layer, followed by the AI2O3 
inter-layer dielectric deposition, (b) Top graphene layer isolation on a separate substrate, followed by transfer onto the bottom 
layer, (c) Top layer Hall bar realization by etching, lithography, metal deposition, and lift-off. The yellow (inner) and red 
(outer) dashed contours in the optical micrograph represent the top and bottom layers, respectively. The scale bars in all panels 
are 10 fim. 




FIG. 2: (color online) (a) Layer resistivities, and (b) densities vs. Vbg measured at T—4.2 K in sample 1. Depending on 
Vbg, both electrons or holes can be induced in the bottom layer; the top layer contains electrons in the available Vbg window, 
owing to unintentional doping. The symbols in panels (a) and (b) represent experimental data, while the lines represent the 
calculated values according to Eqs. (1) and (2). (c) Band diagram across the graphene bilayer heterostructure at Vbg=0 V, 
and (d) at a positive Vbg- Both layers are assumed to be at the charge neutrality point, and aligned with the back-gate Fermi 
level at Vbg = 0. The layers are held at the ground potential, and their thicknesses exaggerated to show the Dirac cones. The 
applied Vbg induces voltage drops Vsio 2 , an d Vai 2 o 3 across the bottom, and inter-layer dielectrics respectively. 



electric capacitance per unit area. Similarly, the Fermi 
energy difference between the two layers is responsible 
for the potential drop across the AI2O3 inter- layer di- 
electric: E F (n B ) = e 2 n T /CAi 2 o 3 + E F (n T ) (2); Cai 2 o 3 
is the AI2O3 dielectric capacitance per unit area. The 
finite Fermi energy of the bottom layer, Ef{ub) in Eq. 
(2), plays the same role with respect to the top layer, as 
the applied Vbg m Eq. (1) for the bottom layer. Equa- 
tions (1) and (2) allow us to determine ub and hit as a 
function of Vbg- This model can be adjusted to include 
finite layer densities at Vbg=0 V. The layer resistivity 
dependence on Vbg can be understood using a Drude 
model pt,b = {nj~ b € -I j 'T,b)~ 1 where /Mr, and ^b are 
the top and bottom layer mobilities, and the layer den- 



sities tit b = y n T b + n 0T os are adjusted to allow for 
finite carrier densities (nor, it-ob) at the charge neutrality 
point. The data of Figs. 2(a) and 2(b) show a good agree- 
ment between the measured layer resistivities and den- 
sities (symbols), and the calculations (solid lines). The 
layer mobilities, determined from Hall measurements, are 
Ms=5,400 cm 2 /Vs, and /i T =4,500 cm 2 /Vs at T=4.2 K. 

Key insight into the physics of the graphene bilayer sys- 
tem can be gained from Coulomb drag measurements^. 
A current (I Drive) flown in one (drive) layer leads to a 
momentum transfer between the two layers, thanks to 
the inter-layer electron-electron interaction. To counter 
this momentum transfer, a longitudinal voltage (Vjj rag ) 
builds up in the opposite (drag) layer. The polarity of 
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FIG. 4: (color online) porag vs. Vbg measured in sample 2 
for T < 77 K. As T is reduced, pDrag exhibits mesoscopic 
fluctuations with increasing amplitude, which fully obscure 
the average drag at the lowest T. The traces are shifted for 
clarity; the horizontal dashed lines indicate Q for each trace. 



FIG. 3: (color online) Coulomb drag in graphene. (a) Layer 
resistivities (pt.b) and porag vs. layer densities (ut,b) for 
sample 2, measured by sweeping Vbg at T=250 K. The bilayer 
probes three different regimes: hole-hole, electron-hole, and 
electron-electron, (b) pDrag vs. Vbg at different T values, 
from 250 K to 77 K (solid lines). Inset: maximum pDrag vs. 
T 2 in the electron-hole and electron-electron regimes. The 
different a;-axis, i.e. ns and nr of panel (a) and Vbg of panel 
(b), apply to both panels. 



Vorag depends on the carrier type in the two layers, and 
is opposite (same) polarity as the voltage drop in the 
drive layer when the both layers have the same (oppo- 
site) type of carriers. The drag resistivity is defined as 
PDrag = {W / L)V Dra g / 1 Drive-, where L and W are the 
length and width of the region where drag occurs. pDrag 
vs. Vbg measured at T=250 K in sample 2 is shown 
in Fig. 3(a), along with the layer resistivities, pr and 
Pb- Unlike sample 1 data (Fig. 2), the charge neutral- 
ity (Dirac) points of both layers can be captured in the 
experimentally accessible Vbg window. Consequently, 
depending on the Vbg value, sample 2 can probe three 
different regimes: a hole-hole bilayer, for Vbg < -15 V, 
an electron-hole bilayer for -15 V < Vbg < -2 V, and 
an electron- electron bilayer for Vbg > -2 V. The depen- 
dence of pb and pr on Vbg of Fig. 3 is also in good 
agreement with the model presented in Fig. 2. Consis- 
tent with the above argument, pDrag is positive in the 
electron-hole bilayer regime, negative in the hole-hole or 
electron-electron regime, and changes sign when either 



the top or the bottom layer are at the charge neutrality 
point. Standard consistency checks^ ensured the mea- 
sured drag is not affected by inter-layer leakage current. 

For two closely spaced two-dimensional systems, when 
the ground state of each layer is a Fermi liquid, and 
the inter-layer interaction is a treated as a perturbation, 
the 

PDrag depends on layer density (n) as oc 1/n 3 ' 2 , 
on temperature as oc T 2 , and inter-layer distance (d) 
as oc l/d 4 j2 Likewise, the Coulomb drag resistivity in 
graphene, calculated in the Fermi liquid regime using 
Boltzmann transport formalism and the random phase 
approximation for the dynamic screening is- 8 : PDrag = 

~~tt^H^ — 3/f 3 /2 (3); ks is the Boltzmann con- 

e a n ^ n ^ 

stant, C(3) = 1.2, e and is the dielectric permittivity. A 
separate effect which has been theoretically advanced as 
the representative Coulomb drag mechanism in graphene 
is trigonal warpingiS,. Figure 3(b) data shows pDrag vs. 
Vbg measured for sample 2 for T values between 77 K 
and 250 K. Away from the bottom layer charge neutral- 
ity point, the pDrag magnitude decreases with increasing 
riB and nx- A power law, pDrag oc l/(n B n^) fitting to 
Fig. 3 data for Vbg > yields an exponent a = 1.25 
±0.25, which depends little on temperature. We note 
that the magnitude of pDrag is a factor of ~ 10 2 lower 
than the values expected according to Eq. (3). While 
further theoretical work is needed to explain this dis- 
crepancy, a possible explanation is that Eq. (3) is valid 
for high densities and/or large inter-layer spacing such 
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FIG. 5: pDrag vs. B measured at T—0.3 K in sample 3, 
showing mesoscopic fluctuations similar to Fig. 4 data. Inset: 
pDrag vs. B data autocorrelation. 



that kp • d ^> 1 (&f denotes the Fermi wave-vector)^; 
for Fig. 3 data kp ■ d < 3 at all layer densities. The 
PDrag oc (ksT) 2 dependence, which stems from the al- 
lowed phase space where electron-electron scattering oc- 
curs, is followed closely for temperatures between 70 K 
and 200 K (Fig. 3(b) inset), and softens for T > 200 
K. Figure 3 data shows a smooth crossover for pDrag 
through f2, from the electron- hole to the electron- 
electron regime [blue (shaded) corridor of Fig. 3(b)]. The 
crossover can be explained by the co-existence of electron 
and hole puddles near the charge neutrality point of the 
bottom layer, which generate drag electric fields of oppo- 
site sign, and cancel the pDrag- 

A remarkable transition in the drag resistance is ob- 
served for T lower than 50 K (Fig. 4). As T is re- 
duced, the pDrag data starts to develop fluctuations su- 
perposed on the average pDrag vs. Vbg dependence of 
Fig. 3 and Eq.(3), valid for diffusive transport. The 
PDrag fluctuations, which are reproducible in different 
measurements, grow in amplitude as T is reduced, and 
fully obscure the average diffusive drag below 20 K. This 
manifestation of mesoscopic physics at elevated temper- 
atures is a consequence of the phase coherence length 
(L v ) increasing with reducing T, and represents the 
counterpart of universal conductance fluctuations^ in 



Coulomb dra g 21 ' 22 . Figure 4 data reveal that PDrag fluc- 
tuation amplitude reaches a maximum near the charge 
neutrality point of the bottom layer (Vbg = — 1 V), 
and increases, albeit slowly as T is decreased. Theo- 
retical arguments^ 2 - indicate that the drag conductivity 

(cTDrag = PDrag/ PTPb) fluctuation amplitude (So Drag) 

depends on relevant length scales and temperature as 
fi&Drag oc T ■ (L^l) / L; I is the electron mean free path. 
For the temperature range examined in Fig. 4, I can 
be considered constant, as the mobility is weakly depen- 
dent on T. Assuming the electron-electron interaction 
is the main phase-breaking mechanism in graphene 2 ^, 
hence L v = l\J Ep /2ksT, the T dependence of Soorag 
and Sporag should follow a oc T -1 / 2 dependence, in good 
agreement with Fig. 4 data. 

To probe the signature of weak localization in Coulomb 
drag, in Fig. 5 we show an example of pDrag vs. per- 
pendicular magnetic field (B) data, measured in sam- 
ple 3 at T=0.3 K and Vbg—0 V; both layers contain 
electrons with layer densities rip = 1.4 x 10 11 cm" 2 , 
and riB — 1-5 x 10 11 cm~ 2 . Similar to the Vbg de- 
pendence of Fig. 4, the pDrag vs. B data shows re- 
producible mesoscopic fluctuations. The auto-correlation 
function (C(AB)) of Fig. 5 data reveals a correlation 
field B c — 47 mT, which corresponds to a phase co- 
herence length Lp = \fhfeB c = 300 nm. Similar L v 
values have been extracted from ensemble average mea- 
surements using scanning gate microscopy 2 ^. 

In summary, we demonstrate independently contacted 
graphene bilayers, and probe the Coulomb drag in this 
system. At elevated temperatures the drag resistance de- 
pendence on density and temperature are consistent with 
the Fermi liquid theory. At reduced temperatures, the 
drag exhibits mesoscopic fluctuations which obscure the 
average drag, a result of the interplay between electron- 
electron interaction and phase coherent transport. 
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